The intriguing life of massive galaxies 
Proceedings IAU Symposium No. 295, 2013 
D. Thomas, A. Pasguali & I. Ferreras, eds. 



©2013 International Astronomical Union 
DOI: 00.0000/X000000000000000X 



Assembly Histories and Observational 
Properties of Simulated Early-type Galaxies 

Peter H. Johansson 1 

1 Department of Physics, University of Helsinki, 
Gustaf Hallstromin katu 2a, FI-00014 Helsinki, Finland 
email: Peter . Johansson@helsinki . f i 

Abstract. We demonstrate that massive simulated galaxies assemble in two phases, with the 
initial growth dominated by compact in situ star formation, whereas the late growth is dominated 
by accretion of old stars formed in subunits outside the main galaxy. We also show that 1) 
gravitational feedback strongly suppresses late star formation in massive galaxies contributing to 
the observed galaxy colour bimodality that 2) the observed galaxy downsizing can be explained 
naturally in the two-phased model and finally that 3) the details of the assembly histories of 
massive galaxies are directly connected to their observed kinematic properties. 
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1. The two phased formation of early-type galaxies 

High-resolution numerical zoom-in simulations of massive early-type galaxies have 
shown that there are two distinct phases in their formation histories (Naab et al. 2007; 
Oser et al. 2010; Lackner et al. 2012). At high redshifts of z ~ 3 — 6 the galaxies assemble 
rapidly through compact (r < r e s) in situ star formation. The later growth of the galax- 
ies proceeds predominantly through the accretion of stars formed in subunits outside the 
main galaxy. The majority of the accreted stars are added to the outskirts of the galaxies 
at larger radii (r > r c ff ) providing a satisfactory explanation for the observed size growth 
of massive early-type galaxies (ETGs) since z ~ 3 until the present-day (Naab et al. 
2009; Bezanson et al. 2009; Oser et al. 2012). 

Here we show that in addition to the size growth of ETGs, the two phased formation 
picture can also shed light on other observational results of ETGs, such as the observed 
galaxy bimodality, the downsizing of massive galaxies and the observed dichotomy of the 
kinematic properties of ETGs. The simulation sample we use to address these questions 
contains 9 galaxies simulated at high-resolution (Johansson et al. 2012) using the Gadget- 
2 code (Springel 2005). The simulations include cooling for a primordial composition, star 
formation and feedback from type II supernovae, but exclude supernova driven winds and 
AGN feedback. We run three models (A2,C2,E2) at very high spatial (e+ = 0.125 kpc) 
and mass resolution (to* ~ 1O 5 M ) with the remaining six simulations (U,Y,T,Q,M,L) 
simulated at a somewhat lower resolution of (e* = 0.25 kpc, to* ~ 10 6 Af Q ). We or- 
ganize the galaxies in our simulation sample into three groups of three galaxies each 
depending on whether their late assembly history (z < 2) is dominated primarily by 
dissipationless minor merging (mostly accreted stars: galaxies C2,U,Y), a mixed dissi- 
pationless/dissipational (mostly accreted and some in situ stars: galaxies A2,Q,T) or a 
primarily dissipational formation history (significant amount of in situ stars: galaxies 
E2,M,L). We find that this classification of the assembly histories of ETGs is useful as it 
broadly separates more massive ETGs forming dissipationlessly from less massive ETGs 
with a more dissipational formation history. 
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Figure 1. The net heating (solid line) and net cooling rates (dashed lines) for virial (r < r v - 1T ) 
non-starforming diffuse gas, for which the density is below the star formation threshold of 
nth < 0.205 cm -3 shown at redshifts of z — (left) and z — 3 (right). The fraction / s f of dense 
starforming gas (n > nth) is also given. Typically, the heating rate dominates over the cooling 
rate at all redshifts for the low-density non-starforming gas. 

2. The bimodality of the local galaxy population 

Recent large-field sky surveys have unequivocally demonstrated that local galaxies 
show a bimodal distribution. Local galaxies with stellar masses above a critical mass of 
M cr it — 3 x 10 10 M Q are typically red spheroidal galaxies with old stellar populations, 
whereas galaxies below this critical mass are typically blue, star-forming galaxies with 
somewhat younger stellar populations. The observed bimodality is usually explained 
theoretically using models in which the star formation is efficiently quenched in massive 
haloes above M ~ 10 12 M Q (Dekel & Birnboim 2006; Gabor et al. 2011). 

The late assembly of our simulated galaxies is dominated by dry minor merging build- 
ing up the accreted stellar component. The infalling satellite galaxies captured through 
dynamical friction will cause gaseous wakes from which energy is transferred to the sur- 
rounding gas. Collectively together with the heating caused by supersonic collisions and 
shocks caused by infalling cold gas this process is deemed gravitational heating. In Fig. 
[T] we compare the shock-induced heating rates of the diffuse non-starforming gas with 
the corresponding cooling rates. The heating rates dominate over the cooling rates at all 
redshifts with the more massive galaxies (U,A2) showing higher heating rates than the 
slightly lower mass galaxies (C2,E2), as expected by the scaling of the gravitational feed- 
back energy, (Ai?) grav oc v^, where v c is the circular velocity of the galaxy (Johansson 
et al. 2009b) . The inclusion of gravitational heating helps in maintaining a hot gaseous 
halo and thus inhibits star formation contributing to the observed galaxy bimodality. 
However, gravitational heating is predominantly important in the outer parts of galaxies 
and some form of additional feedback, most probably AGN feedback (e.g. Johansson et 
al. 2009a), is required to stop late central star formation in very massive galaxies (e.g. 
galaxy U at z = 0). 



3. Downsizing of massive galaxies 

Several recent observations have shown that old, massive red metal-rich galaxies were 
already in place at high redshifts of z ~ 2 — 3. This observational result can be seen as 
a manifestation of cosmic downsizing, in which galaxies seem to form anti-hierarchically 
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Figure 2. The stellar mass assembly histories of our simulated galaxy sample, with the solid 
colour bars showing the contribution of in situ formed stellar mass and the dashed colour bars 
representing the contribution of accreted stellar mass shown at redshifts of z = (left) and z — 3 
(right). At high redshifts (z > 3) the galaxies assemble rapidly through in situ star formation, 
whereas the late (z < 3) assembly history is dominated by accreted stars, with the more massive 
galaxies ending up with a proportionally larger fraction of accreted stars. 



in the sense that the most massive galaxies formed a significant proportion of their stars 
at high redshifts, compared to lower mass systems that exhibit a more continuous star 
formation history throughout the cosmic epoch (e.g. Glazebrook et al. 2004). 

In Fig. [2] we show the assembly histories of all our simulated galaxies depicting sep- 
arately the masses in the in-situ formed and accreted stellar components. At z = 3 the 
stellar components in all galaxies have been assembled rapidly through mainly in situ 
star formation, fueled by cold gas flows and hierarchical mergers of multiple star-bursting 
subunits. At lower redshifts (z < 3) the subsequent growth of the stellar component pro- 
ceeds predominantly through the accretion of existing stellar clumps. The galaxies in Fig. 
[5] are ordered in decreasing final stellar mass from top to bottom and we can immediately 
see that the fraction of accreted stars at z — increases as a function of galaxy mass. 
The accreted stellar component forms on average in low mass galaxies very early in the 
Universe. However, the accreted stellar component is added to the massive galaxies much 
later and at lower redshifts than the in situ formed stars. In addition, the metallicity of 
the central in situ formed stars is on average higher than the metallicity of the accreted 
stars that are formed in lower mass galaxies and added later to the outskirts of galaxies, 
resulting in a negative metallicity gradient, in agreement with the observations. Thus, 
the counter-intuitive concept of downsizing can be explained in the two phased forma- 
tion mechanism. Massive galaxies form their central stellar mass in situ and then accrete 
substantial amounts of stars that were formed even earlier in smaller subsystems. Hence 
by z ~ 2 — 3 the most massive galaxies have the oldest stellar populations compared to 
lower mass galaxies that are still forming in situ stars. 



4. Kinematic properties of massive galaxies 

Observations in the late 1980s using slit spectroscopy found that the ETG popu- 
lation can be broadly separated into a population of massive slowly-rotating systems 
(v/a < 0.1) with boxy isophotes and a population of fast-rotating (v/a ~ 1) ETGs with 
more disky isophotes found typically at somewhat lower masses. Recent results from the 
volume limited ATLAS 3D survey utilizing a modern integral-field-unit (IFU) confirmed 
this dichotomy showing that about ~ 15% of local ETGs rotate slowly with no indica- 
tions of an embedded disk component, whereas the majority (~ 85%) of the local ETGs 
show significant disk-like rotation (Cappellari et al. 2011). 




Figure 3. The effective ellipticity (e e ff) is plotted against the ratio of the major axis rotation 
and central velocity dispersion (« roa j / Co) ■ The contours show the 95% probability location of the 
galaxies in the (e e fi — v m ajA r o)- plane derived from 500 random viewing angles of our simulation 
data. The overplotted symbols are observational data from (Bender et al. 1994) demonstrating 
that our simulated galaxies seem to be largely consistent with the observations. The dashed line 
shows the theoretical value for an oblate isotropic rotator. 

We derive the kinematic properties of our simulated galaxies using 500 random pro- 
jections in order to assess the mean properties of the simulated galaxies averaged over 
all sightlines. In Fig. [3] we plot the 95% probability of finding a simulated galaxy in the 
£efi — (fmaj/co) plane, where e e ff measures the effective ellipticity of the galaxies and 
the ratio of the major axis rotation and velocity dispersion (v m aj/co) is a measure of 
the rotational support of the galaxies. Galaxies U,C,E2 show projections with very low 
rotational support, whereas the other galaxies (Y,A2,Q,T,L,M) show significantly more 
rotational support in very good agreement with the theoretical prediction for an oblate 
isotropic rotator shown by the dashed line in Fig. [3l The most massive galaxy U shows 
mostly round projections with very low e c g, whereas all the other galaxies show projec- 
tions extending in e e ff from zero up to 0.4. We see a correlation between the rotational 
support and the in situ/accreted fraction, both galaxies U,C that have a large accreted 
fraction are slowly rotating and galaxy Y would most probably also been slowly rotating 
if it had not experienced a late (z < 0.5) major merger. Thus, within our rather narrow 
mass range the majority of the simulated galaxies are consistent with being rotationally 
supported disky ellipticals, whereas the most massive galaxy in our sample is consistent 
with being a roundish slow-rotator. 
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